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Abstract 

Motivated by recent hints in particle physics and cosmology, we study the realization of eV-scale 
sterile neutrinos within both the seesaw mechanism and flavor symmetry theories. We show that 
light sterile neutrinos can rather easily be accommodated in the popular ^4 flavor symmetry 
models. The exact tri-bimaximal mixing pattern is perturbed due to active-sterile mixing, which 
we discuss in detail for one example. In addition, we find an interesting extension of the type I 
seesaw, which can provide a natural origin for eV-scale sterile neutrinos as well as visible admixtures 
between sterile and active neutrinos. We also show that the presence of sterile neutrinos would 
significantly change the observables in neutrino experiments, specifically the oscillation probabilities 
in short-baseline experiments and the effective mass in neutrino-less double beta decay. The latter 
can prove particularly helpful in strengthening the case for eV-scale sterile neutrinos. 
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I. INTRODUCTION 



Neutrino oscillation experiments have provided firm evidence that, in contrast to the 
prediction of the Standard Model (SM), neutrinos are massive and that their flavors change 
during propagation. Their unusual mixing pattern and the smallness of their masses makes 
the explanation of the origin of neutrino masses and leptonic flavor mixing one of the most 
challenging problems in particle physics. In the standard neutrino oscillation picture three 
active neutrinos are involved, with mass-squared differences of order 10~^ and 10"'^ eV^. 
The absolute mass scale of neutrinos is also constrained to be less than around 1 eV from 
tritium beta decay experiments as well as cosmological observations. For a recent review on 
our current understanding of neutrino observables, see Ref. [1]. 

Despite the successful achievements of solar, atmospheric, reactor and accelerator neu- 
trino experiments, there are experimental anomalies that cannot be explained within the 
standard three neutrino framework. In particular, the possible presence of sterile neutrinos 
points towards non-standard neutrino physics. The issue of the LSND and MiniBooNE 
results has been around for some time, and is frequently interpreted as a hint towards the 
presence of one or two sterile neutrino states [2-5] . Recently this debate has been re- ignited 
by a reevaluation of the anti-neutrino spectra of nuclear reactors [6], which leads to increased 
fluxes. As a result, the negative results of previous reactor experiments can in fact be inter- 
preted as the observation of a flux deficit, and this in turn can be explained by additional 
sterile neutrinos with masses at the eV scale [6, 7]. Interestingly, current results from pre- 
cision cosmology and Big Bang nucleosynthesis mildly favor extra radiation in the Universe 
beyond photons and ordinary neutrinos. While this could be any relativistic degree of free- 
dom, the interpretation in terms of additional sterile neutrino species is straightforward. 
Indeed, several cosmological parameter fits (e.g., the analyses of the CMB or SDSS data 
sets in Ref. [8]) have been found to be compatible with more radiation than that predicted 
by the SM particle content. This is supported by the recently reported higher '^He abun- 
dance [9, 10], which in the framework of Big Bang nucleosynthesis can be accommodated 
by additional relativistic degrees of freedom, as it leads to earlier freeze-out of the weak 
reactions, resulting in a higher neutron-to-proton ratio. It is rather intriguing that hints 
for the presence of sterile neutrinos are given by fundamentally different probes: neutrino 
oscillations. Big Bang nucleosynthesis and observation of the Universe's structure. This is 
the situation which motivates the present study. 

Sterile neutrinos, if they exist, would lead to rich experimental phenomena. By definition, 
sterile neutrinos do not directly enter the weak interactions. However, their admixture with 
active neutrinos would modify the neutrino flavor mixing and lead to observable effects in 
neutrino oscillation experiments. Furthermore, due to that admixture they could interact 
with gauge bosons, resulting in significant corrections to certain non-oscillation processes, 
e.g., in the neutrino-less double decay (Oz//3/3) amplitude [11-13] or in beta decay spectra, 
such as in the KATRIN experiment [14, 15]. It has also been pointed out that an eV-scale 
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sterile neutrino would significantly affect the the atmospheric neutrino fluxes in the energy 
range 500 GeV to a few TeV, and these effects could be studied in the IceCube detector [16- 
18]. Other aspects that are in principle modified by such sterile states are supernova physics 
[19], solar neutrinos [20], or the interpretation of cosmological data [21]. We will add to 
this discussion by updating the predictions for Ou/J/S with the recent results for the relevant 
mixing parameters. We consider all possible neutrino mass spectra, namely four in case of 
one sterile neutrino, and eight in case of two sterile neutrinos. For the latter case we also 
point out some properties of short-baseline oscillation probabilities in four often overlooked 
possible spectra. 

On the other hand, we consider the non-trivial theoretical origin of eV-scale sterile neutri- 
nos. In what regards neutrino masses and their peculiar mixing structure, the key words are 
seesaw mechanism [22-25] and flavor symmetries (see [26, 27] for recent reviews). We will 
show that popular models based on the A4 flavor symmetry can easily be modified to take 
sterile neutrinos into account. As a result, the mixing structure of the original model, tri- 
bimaximal mixing (TBM) in our example, is perturbed, and corrections to the tri-bimaximal 
values arise. This will be a general feature of such approaches, but here we focus on one 
concrete example, which we will outline in detail.^ 

Sterile neutrinos are a necessary ingredient of the canonical type I seesaw mechanism, 
though they are naturally assumed to be many orders of magnitude heavier than the SM 
scale of 10^ GeV. As a result, their mixing with the SM particles is highly suppressed. If 
one brings one of the heavy states down to the eV scale, it is possible to generate a sterile 
neutrino with the correct mass and (potentially) the correct mixing with the SM leptons 
to explain the data. However, another more non-trivial case can also be studied, namely 
an extension of the type I seesaw with additional heavy (i.e. heavier than the SM scale) 
neutral fermions. This approach generates sterile neutrinos without the need to initially 
have states with eV masses, and is more in the seesaw spirit. 

The remaining parts of this work are organized as follows: in Sec. II, we outline the for- 
malism of neutrino mixing in the presence of sterile states and summarize the phenomeno- 
logical consequences of sterile neutrinos in neutrino-less double beta decay. Short comments 
on short-baseline neutrino oscillations are delegated to the Appendix. Then, in Sec. Ill, 
we discuss how to embed sterile neutrinos into one particular A4 flavor symmetry model, 
and describe the resulting deviations from exact TBM. Section IV is devoted to a general 
overview of the realization of sterile neutrinos in seesaw models. Finally, we conclude in 
Sec. V. The individual sections are largely independent of each other, but their results and 
methods could be combined. However, we feel that a separate discussion of each aspect is 
more suitable for the present discourse. 



Note that a general analysis of the mixing of active and sterile neutrinos has been presented in Ref. [28]. 
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II. PHENOMENOLOGICAL CONSEQUENCES OF STERILE NEUTRINOS 



In this section we will outline some properties of neutrino parameters and evaluate the 
contributions to the effective mass relevant for neutrino-less double beta decay. 

A. Neutrino mixing with sterile neutrinos 

In the presence of = n — 3 sterile neutrinos, the neutrino mass matrix is an n x n matrix 
m,y, which can be diagonalized by means of an n x n unitary matrix U. The neutrino flavor 
eigenstates uj (for / = e, /i, r, Si, S2, ■ ■ ■ , Sn-3) are then related to their mass eigenstates Ui 
(for i = 1,2, . . . ,n) via 



(1) 



In general, for n massive families including = n — 3 7^ massive sterile neutrinos, one has 
n — 1 = Ug + I Majorana phases, 3(n — 2) = 3(72^ + 1) mixing angles and 2n — 5 = 272^ + 1 Dirac 



phases. The number of angles and Dirac phases is less than the naive \ 



n[n 



1) angles and 



\{n — l){n—2) phases, because the \ 



nJn^ 



■1) rotations between sterile states are unphysical. 



For illustration, in the case of only one sterile neutrino, U is typically parameterized by 

U = -R34-R24-R14-R23-R13-R12-P , (2) 
where the matrices Rij are rotations in ij space, i.e. 
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(3) 



where Sij = sin 6ij, = cos%. The diagonal P matrix contains the three Majorana phases 
a, (5 and 7: 

P = diag (1, e''^'\ e'(/5/2+'5"), e*(V2+<5i4)) _ (4) 

Note that there are in total three Dirac CP- violating phases The above definition of 
P is constructed in such a way that only Majorana phases show up in the effective mass 
governing neutrino-less double beta decay (see below). Similarly, one can parameterize the 
mixing matrix for 2 sterile neutrinos as 



U — R25R34R25R24R23R15R14R13R12P 



(5) 
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where P = diag(l, e*"/^ e'^^/^+^i^)^ ei(7/2+5i4)^ ei(0/2+5i5))_ 

The mass-squared differences associated with the current hints for sterile neutrinos are 
much larger than the ones responsible for solar and atmospheric oscillations. The sterile 
neutrinos are thus separated in mass from the active ones and can either be heavier or 
lighter than the active ones. The active neutrinos can be normally {m^ > m2 > mi) or 
inversely (m2 > rrii > m^) ordered. Therefore, if there is one sterile neutrino, there are in 
total four possible mass orderings. Our nomenclature for the schemes is such that we go 
from top to bottom, i.e. if ^ "t-i 2,3, we denote this 1+3 scenario as "SN" if the active 
neutrinos are normally ordered and "SI" if they are inversely ordered. Analogously, if the 
sterile state is lighter than the active ones (r«i,2,3 ^ i^s, 3+1 scenario), we name the schemes 
"NS" and "IS", respectively. Obviously, 3+1 scenarios are less attractive since they induce 
(more) tension with the cosmological bound on the sum of neutrino masses. The models we 
will present in later sections indeed predict 1+3 scenarios. 

In the case of two sterile neutrinos with masses m^^ and m^j, there are three classes 
of mass spectra: ms^^nis^ ^ "21^2,3 (2+3 scenarios), mi^2,3 ^ ms^.ms^ (3+2) and [13] 
^ ?Tii^2,3 ^ (1+3+1). The latter class of 1+3+1 orderings, first noted in Ref. [13], 
has received less attention in the past, and was for instance not taken into account in the 
fits of Refs. [2-4]. The 1+3+1 cases have the interesting property that the active states 
are sandwiched between the two sterile ones. Current oscillation data cannot distinguish 
1+3 from 3+1, or 3+2 from 2+3 scenarios, but are sensitive to 1+3+1 vs. 2+3/3+2 (see 
the Appendix). Ref. [7] fits the world's short-baseline data also within the 1+3+1 scenario 
and finds that the fit is slightly better than that for 2+3/3+2. The 2+3 cases are denoted 
"SSN" or "SSI", while the 3+2 cases are called "NSS" or "ISS". In what regards the 1+3+1 
cases, two possible permutations of the mass spectrum should be distinguished. One of the 
sterile neutrinos is separated from the active ones by a larger mass gap, associated with 
the larger of the two mass-squared differences. If the heavier of the two sterile neutrinos is 
connected with the larger mass-squared difference, the scenarios are denoted "SNSa" and 
"SISa" , respectively; if the heavier sterile state is connected with the smaller mass-squared 
difference, we call the cases " SNSb" and "SISb" . The 2+3 scenarios are more attractive than 
the others because they predict a smaller sum of masses. The individual masses mi^2,3,4,5 
expressed in terms of the mass-squared differences Amg, Am^, Am\i and Amg^ can be 
found in Ref. [13], where the generalization to three sterile neutrinos (with sixteen possible 
mass orderings) is also discussed. 

Table I shows the best-fit and 2a ranges of the relevant parameters used for the analysis 
in this work. The best-fit values are taken from the global fit in Table II of Ref. [7]. In 
their analysis of the 3+1/1+3 scenarios, the authors of Ref. [7] find several different allowed 
regions in the ^m\i — sin^^i4 parameter space, at 2a. For convenience we use the region 
around the best-fit point, and (since the ranges are not available) allow the parameters in 
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TABLE I: Best-fit (from Ref. [7]) and estimated 2a values of the sterile neutrino parameters. 





parameter Am4]^ [eV] 


|t/e4p 


01 


[eV] 




3+1/H 


best-fit 


1.78 


0.023 








1o 


1.61-2.01 


0.006-0.040 








3+2/24 


best-fit 


0.47 


0.016 


0.^ 


^7 


0.019 




0.42-0.52 


0.004-0.029 


0.77- 


-0.97 


0.005-0.033 


1+3+1 


best-fit 


0.47 


0.017 


0.^ 


^7 


0.020 


2a 


0.42-0.52 


0.004-0.029 


0.77- 


0.97 


0.005-0.035 



the 3+2, 2+3 and 1+3+1 cases to have 2a uncertainties of the same relative magnitude as 
those in the 3+1/1+3 scenario. The data favor the presence of two sterile neutrinos, mostly 
because they allow different neutrino and anti-neutrino probabilities, thus alleviating the 
tension between the LSND and MiniBooNE results. As mentioned above, 1+3+1 scenarios 
have a slightly better fit than 3+2/2+3 cases. 

We note in addition that the recent results from the T2K [29] and MINOS [30] experiments 
strengthen the existing hints [31-33] for non-zero ^13, and that the analysis in Ref. [34] finds 
evidence for ^13 > at the level of > 3(t. It is not yet evident whether these new data will 
improve or worsen the fits in the various sterile neutrino scenarios; on the other hand it can 
also be argued that the T2K result is not due to 6*13 but is actually another signature of 
sterile neutrinos [35]. 



B. Neutrino-less double beta decay 

Neutrino-less double beta decay (Oi^/3/3) is the only realistic test of lepton number vio- 
lation. While there are several mechanisms to mediate the process (e.g. heavy neutrinos, 
right-handed currents or SUSY particles), light Majorana neutrino exchange is presumably 
the best motivated scenario. We will work in this standard interpretation of OufSP, and 
study the effect of massive sterile neutrinos in the eV range. This section is an update of 
the results in Refs. [12, 13]; we also use the best-fit values in Table I to report approximate 



numerical lower limits for the complementary mass observables = y |[/ej|^m? and "^rrii, 
constrained by beta decay and cosmology, respectively. 

In the presence of one sterile neutrino, the effective neutrino mass in Ouf^P is given by 

/ \ 1222 ,222 iai22 iS , 2 i'yl //?\ 

("^ee)4^. = \c12Cr3Cum1 + Si2Ci3Ci4m2e + Sy^c^^iTise " + s^^mie^\ , (6) 

using the parameterization in Eq. (2). If the sterile neutrino is heavier than the active ones, 
the approximation 
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1+3, Normal, SN 



1+3, Inverted, SI 




3 V (best-fit) 
3 V (2a) 
1+3 V (best-fit) 
1-1-3 V (2a) 



0.001 



0.01 




m,- „ (eV) 

light ^ ' 



3+1, Normal, NS 



3+1, Inverted, IS 




FIG. 1: The allowed ranges in the {Triee) ~ ^'Ught parameter space, both in the standard three- 
neutrino picture (unshaded regions) and with one sterile neutrino (shaded regions), for the 1+3 
(top) and 3+1 (bottom) cases. 



holds, where {mf.^^y = c\2c\-^^i + s\2c\2jn2e^" + s\^m^e^l^ is the standard expression for three 
active neutrinos. The upper panel of Fig. 1 displays the allowed range of {mee){i+-i)u as a 
function of the lightest mass mught, using data from Refs. [7, 33]. Also shown in this plot 
and the following ones is the allowed range of \{mee)zu\i i-e. the standard case. Its crucial 
features (see [36] and references therein) are that |(mee)3;/| can vanish exactly in the normal 
hierarchy case, but cannot vanish in neither the inverted hierarchy nor the quasi- degenerate 
case. This standard behavior can be completely mixed up by the presence of one or more 
sterile neutrinos: if the lightest neutrino mass is reasonably small, e.g. mnght < 0.01 eV, the 
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allowed range of {mee){i+3)u is dominated by the term s^4y^Am|^ ^ 0.031 eV, which means 
that {friee) {i+3)u cannot vanish in the normal ordering case (the contribution of the two light 
active neutrinos cannot cancel that of the sterile neutrino). However, in the inverted ordering 
("^ee)(i+3)i/ can vanish even for very small active neutrino masses. The effective mass can 
also be zero in the regime where the active neutrinos are quasi-degenerate (mught > 0.1 eV). 

This modified behavior of (mge) is of particular interest in the inverted ordering case, 
since the usual lower bound on the effective mass (cf. the solid and dashed lines in Fig. 1) is 
no longer valid. If future Oz//3/3 experiments measure a tiny effective mass and the neutrino 
mass hierarchy is confirmed to be inverted from long-baseline neutrino oscillations, the sterile 
neutrino hypothesis would be an attractive explanation for this inconsistency. In the 3-1-1 
case, and ^ rrii are dominated by the sterile contribution, i.e. ^ ^/\Ue^\^Am^ ~ 
0.2 eV and X]"^* ~ V ^"^li — 1-3 eV, respectively. 

The lower panel of Fig. 1 shows the effective mass when the sterile neutrino is lighter than 
the active ones (the 3-1-1 scenario). In this case there are three quasi-degenerate neutrinos at 
the eV scale, with their mass given by ^^/Am^ ~ 1.3 eV (this mass scale governs predictions 
for m^). The effective mass then takes its standard form for quasi-degenerate neutrinos: 



(mee)(3+i)!. ^ Y Am|i Y 1 - sin 26^12 sin a/2 . (8) 

However, this situation is relatively disfavored by cosmological bounds on the sum of 
neutrino masses [8], since ^ ~ 4 eV. If taken at face value, the current 

limit of about 0.5 eV on the effective mass means that v^l - sin^2^i2 sin^a/2 < 0.4, thus 
already putting strong constraints on the solar neutrino mixing angle and in particular the 
Majorana phase. 

If there are two sterile neutrinos, the effective mass reads 

/ \ 12222 ,2222 joi222 iB , 2 2 i'y , 2 id)\ 

(9) 

with the additional Majorana phase. In the 2-1-3 cases where both of the sterile neutrinos 
are at the eV scale, (mee) can be approximated by 



(2+?,)u 



^51 



(10) 



in analogy to the 1-1-3 case. The upper panel of Fig. 2 shows the allowed regions in this 
case: the phenomenology is similar to that discussed for l-|-3 above, except that the pres- 
ence of two sterile terms in Eq. (10) allows (mge) to take smaller values in the hierarchi- 
cal region for the normal ordering. The inverted ordering is essentially the same as in 
the 1-1-3 case. The two eV-scale sterile neutrinos dominate the KATRIN observable, with 
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2+3, Normal, SSN 2+3, Inverted, SSI 




3+2, Normal, NSS 3+2, Inverted, ISS 




FIG. 2: Same as Fig. 1, for the 2+3 (top) and 3+2 (bottom) cases. 



^ A/|t/e4P ^"^41 + l^^esP^^^si — 0.16 eV, and the sum of masses is approximately 
Y^mi> y^Am^ + ^/Am^ ^ 1.6 eV. 

The lower panel of Fig. 2 displays the 3+2 cases, where the sterile neutrinos are lighter 
than the active ones. Since the three active neutrinos are quasi-degenerate at the eV scale, 
with their mass given by the largest sterile mass-squared difference ^y Aml^ ~ 0.69 eV, 
Eq. (8) applies for {mee)(3+2)u and the value of is also set by this scale. The mass 
ordering of the active states plays no role. Once again these scenarios are disfavored by 
cosmology, since ^ mj ^ 3^/Am^ + y^Am|7^^'Am|^ ~ 3.4 eV. 

For completeness we include plots of the effective mass in the 1+3+1 schemes, where 
the active neutrinos are sandwiched between the sterile ones. The plots in Fig. 3 show 
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1+3+1, Normal, SNSa 1+3+1, Inverted, SISa 




m, „ (eV) 

light ^ ' 

FIG. 3: Same as Fig. 1, for the 1+3+la (top) and 1+3+lb (bottom) cases. 

that there is httle difference between these scenarios and the 3+2 cases. The main differ- 
ence arises in SNSb and SISb scenarios, in which the effective mass is governed by three 
quasi-degenerate neutrinos with mass -^/Arn^ ^ 0.93 eV, plus a (coherent) contribution of 
|[/e4p \/ Am4]^ + Amg]^ ~ 0.03 eV. The KATRIN observable is also dominated by this mass 
scale, and the sum of masses is ^ mj ^ 3y^Am|j" + y^AmJ~+"^"^ — 4.0 eV. 

The results of this section will not be significantly modified in the presence of new 
data [29] on sin^6'i3. In the hierarchical region with normal ordering, the lower limit on the 
effective mass decreases with sin^^is, so that the larger best-fit value from Ref. [34] would 
have a small effect on (mge) in the standard three-neutrino scenario. However, in the sterile 
case the dominant terms in Eqs. (7) and (10) are independent of sin^^i3. 



10 



As discussed in the different cases above, if the predictions of the effective mass are above 
0.2 eV, which is true for most of the possible scenarios, then KATRIN will see a signal as 
well. Indeed, Refs. [14, 15] show that sterile neutrinos with masses and mixings in the 
ranges considered above will be observable in KATRIN, and will be distinguishable from 
active neutrinos. This effect should be seen in almost all cases; the smallest contribution to 
the observable m^? arises in the SSN and SSI cases. 

In what regards the sum of masses, the smallest result holds again for SSN and SSI 
(^mj ^ 1.6 eV); the value of exceeds this in all other cases. The current upper 

bound from cosmology is approximately Y^rrii ^ 0.5 eV [37], which introduces tension with 
the sterile neutrino scenarios presented above. However, this measurement is highly model 
dependent and depends on the cosmological data set used. The Planck satellite should 
in the near future be able to reach a sensitivity of ~ 0.1 eV, thus providing more precise 
limits on the absolute neutrino mass. In addition, the present hints for extra radiation from 
analysis of CMB anisotropy measurements will be further constrained by Planck data, with 
a precision of AN^s ~ ±0.3 [8, 38], where A^efr is the effective number of thermally excited 
neutrino degrees of freedom. This means that Planck could in fact rule out a sterile neutrino 
species before beta decay or 0^(3/3 experiments see a signal. 

III. STERILE NEUTRINOS IN FLAVOR SYMMETRY MODELS: AN EXAM- 
PLE 

In this section we modify a popular flavor symmetry model, which predicts tri-bimaximal 
mixing and is based on the group yl4, in order to accommodate one or two light sterile 
neutrinos. 

A. One sterile neutrino 

In the Altarelli-Feruglio model from Ref. [39], neutrinos get mass from effective oper- 
ators, and the judicious choice of particle and fiavon content along with the correct vacuum 
expectation value (VEV) alignment leads to TBM. The relevant particle assignments are 
shown in Table II. Note the presence of an additional symmetry to separate the neutrino 
and charged lepton sectors, and the Froggatt-Nielsen f/(l)FN to generate a hierarchy for 
the charged lepton masses. We have also included a sterile neutrino with appropriate 
quantum numbers under the symmetries. 

Leaving Ug aside for the moment, these particle assignments, along with the multipli- 
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TABLE II: Particle assignments of the model, modified from Ref. [39] to include a sterile 
neutrino Ug. The additional Z3 symmetry decouples the charged lepton and neutrino sectors; the 
C/(1)fn charge generates the hierarchy of charged lepton masses and regulates the scale of the 
sterile state. 



Field 
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hu4 ^ ^' 






SU{2)l 
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2 1 1 
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1" 
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1 3 3 
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1 












1 1 u 
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f/(l)F7V 
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1 
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cation rules (see e.g. Refs. [26, 27]), lead to the Lagrangian 



+ ^^{LhuLhu) + ^{(f'LhuLhu) + h.c. + . . . , 



(11) 



where A is the cut-off scale and the dots stand for higher dimensional operators. The 
notation is such that two fields a and b written as {ab) transform as 1, etc. If one chooses 
the real basis for ^4, along with the fiavon VEV alignments^ (C,) = u, (ip) = {v, 0, 0) and 
{if') = {v',v',v'), then the charged lepton mass matrix is diagonal, and the neutrino mass 
matrix 

/n' -L 2d. _^ _^ \ 
/ a -h 3 3 3 

(12) 



is diagonalized by the TBM matrix. 



V 



3 3 3 

M. n' — ^ 

3 " 3 
2d' 
3 



U' 



TBM 



' V3 



\ 



\ L ±. ±. I 

\ Ve V3 y/2 / 



(13) 



In this case My is form-diagonalizable: the eigenvectors are independent of the parameters in 
the neutrino mass matrix. The charged lepton mass hierarchy is generated by the Froggatt- 
Nielsen (FN) mechanism: the [/(1)fn charges 0, 1 and 3 are assigned to the right-handed 
charged singlets r"^, /i*^ and e'^, respectively, and a fiavon Q that carries a negative unit of 
this charge is introduced, suppressing each mass term by powers of the small parameter 



^ Several methods have been proposed to explain the different alignments and solve the so-called "vacuum 
alignment problem" ; these will not be discussed here. Note that we do not introduce more flavons to the 
original model, so that we can assume that the solution of the alignment problem will not be modified. 
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(^)/A = A < 1. Explicitly, 

rria = VaVdj)^^" , (14) 

where Fa is the relevant f/(l)FN charge. 

By assuming that (i) the Yukawa couplings tja, Xa and Xd remain in a perturbative 
regime; (ii) the flavon VEVs are smaller than the cut-off scale and (iii) all ffavon VEVs fall 
in approximately the same range, the authors of Ref. [39] obtain the relation 

v' V u , . 

0.004 <-^-^-<l, 15 
AAA ^ ' 

with the cut-off scale A ranging between 10^^ and 10^^ GeV, and ~ 174 GeV. 

In order to accommodate sterile neutrinos in this framework, we have added an additional 
sterile singlet, Vs-, transforming as 1 under and 1 under Z^. The invariant dimension-5 
operator j^{(p'Lhu)i^s is not allowed by the Z3 symmetry, and we are left with the terms 

X X f 

C-Y, = -^^if' Lhu)us + j^{ip'ip'LK)v, + msV^y, + h.c, (16) 

where is a bare Majorana mass. As we will see, the chosen FN charge forces to be 
at the desired eV scale. The modified 4x4 mass matrix is (cf. Eq. (12)) 



/a + f -i -f e\ 



M' 



where a = 2xa^, d = 2xcr^ and e = y/2xe^^j:^ have dimensions of mass. Note that the 
first three elements of the fourth row of M^^^ are identical because of the VEV alignment 
((/)') = {v',v',v'), which was necessary to generate TBM in the 3 neutrino case; this align- 
ment combined with the /I4 multiplication rules also causes the second term in Eq. (16) 
(proportional to Xj) to vanish. 

The correct mass scales can be obtained by assigning a f/(l)FN charge of F^^ = 6 to Ug, in 
analogy to the mechanism used in the charged lepton sector. Indeed, in order to fit the data, 
the parameters a and d should be between 10~^ and 10~^ eV, the ratio e/m^ ~ (9(10~^) to 
generate sufficient mixing, and ~ 1 eV for the sterile neutrino mass. 

As an explicit example, assume that v'/A ~ v/A ~ u/A ^ 10~^'^, in keeping with the 
constraints from Eq. (15), and that the cut-off scale is A ^ 10^^'^ GeV. In this case one 



4x4 



2d „ _ d 
3 " 3 



2d 

3 ^ 



(17) 



V ■ ■ • msj 



13 



obtains 



e ^ 0.1 



a ~ c? ~ 0.1 



A 



u 



10 



-1.5 



10^1 GeV/ V102 GeV 



J 



10^1 GeV/ V 1011 GeV / V102 GeV 



1012-5 GeV\' 

-1 eV; 



(18) 



A 



with the assumption that the Yukawa couphngs Xa^^e are of order 1, and that A ~ 10~i'^, 
so that {9) is in the same range as the other flavon VEVs. 

The Majorana mass term msU^Ug is doubly suppressed by the ?7(1)fn charge, and there 
are additional terms which can give a contribution to this mass in addition to the bare 
term. From the particle assignments in Table II, the leading and next-to-leading order 
non-vanishing contributions to rris are^ 



V 

A 



u 



3v'^u 
1^ 



(19) 



so that these terms are suppressed by Ai^, and the resulting Majorana mass can be of order 
eV, i.e. 

A V \ 2 ^1012-5 GeV\ 

eV . 



10 



0.5 



10-1.5 



1011 GeV 



) 



A 



(20) 



In this case the contributions from the second and third terms in Eq. (19) are of order 0.1 eV, 
and do not affect the scale of significantly. We conclude that the usual choice of scales 
and charges easily allows a ~ ~ e < in M^^^ from Eq. (17). 

Assuming that the parameters are all real, Eq. (17) is exactly diagonalized by 



U 



2 

V6 


1 K. 

6e N- 





1 ^+\ 

6e N+ 


1 

Vg 


1 

6e 


1 

V2 


1 K+ 
6e N+ 


1 


1 K- 

6e Af_ 


1 
V2 


1 K+ 
6e N+ 





1 
N- 








(21) 



where K-^ 



nis ± \/l2e2 + (a — nisY and A''! 



a—ms:k^l2e^+{a—m,a)'' 



If one 



assumes that a < rris and expands to second order in the small ratio e/rus, the resulting 



^ The term proportional to ^^(ip'ip'ip') vanishes after A4 symmetry breaking. 
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mixing matrix is 





2 


1 





o\ 




/o 








rris 




1 


1 


1 

V2 





+ 











e 

ms 




1 


1 


1 

V2 












rris 





e 


v 











V 




\o 





0/ 



u 



giving the eigenvalues 

mi = a + d , 7712 = a 



+ 



/O 






2m2 
y3e2 









2m2 





\ 






(22) 



3e2 

— , ms = —a + (i . 
ms 



3e2 

7714 = ms H 

m.. 



(23) 



In this case one can see that M^^"^ is not form-diagonahzable anymore: the second and fourth 
column of U are sensitive to the entries of M^^^. The mass-squared differences as well as 
the active-sterile mixing angles sin^ 6*44 (z = 1,2,3) are controlled by the four parameters 
in Eq. (17), with the mixing most sensitive to the ratio e/m^. Note that both the normal 
and inverted orderings are allowed, in contrast to the standard three neutrino version of the 
model, which only allowed the normal ordering. By fitting the parameters a, d, e and m^ 
to the allowed range of the four parameters Amg, Am\, Am^^ and sin^^i4, we find that 
the masses are arbitrary, are not constrained to any particular region, and are in general 
uncorrelated with the mixing parameters. However, the lightest mass increases with sin^^i4, 
which means that the effective mass in OvP/S, (mee) = |a -|- y |, also increases with sin^^i4, 
as expected from Eq. (7). Figure 4 shows the allowed ranges of a,d,e. In general, a and d 
are approximately inversely proportional to each other. 

Comparison of Eqs. (2) and (22), shows that sin ^13 = 0, i.e. this parameter retains its 
TBM value, whereas sin^6'i2 and sin^6'23 receive small corrections: 



sin^^ 



12 





u 


e2 


2 


1 - 




UeA 


2 



1 

^ 3 



21^ 



sin^^. 



If/. 



23 



/i3| 



1- a 



e4| 



u, 



1 1 

2 - 2 



1+ ( — 

m.. 



(24) 



Note the correlation sin^^23 — |(1 ^ sin^6'i2) following from the above expressions. Other 
results of the model are Ugi = Ugs = and f/e4 = U^4^ = 11^4,. The three active-sterile mixing 
angles can be expressed in terms of the matrix elements U ^4 as 



sin^^i4 



\U, 



e4 



sin^^. 



24 







2 


1 - 




u 


e4 


2 



sin^^ 



34 







2 


1 - 




2 _ 




U,A 


2 



(25) 
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a(eV) 

FIG. 4: The allowed values in a — d and a — e parameter space for normal (NO) and inverted (10) 
ordering, obtained by varying each parameter between —0.5 and 0.5 eV, varying mg between —1.5 
and 1.5 eV, and requiring that the oscillation parameters lie in the correct range [7, 33]. 

and the model predicts them all to be of similar magnitude: 

sin^^u ~ sin2^24 ^ sm'^Ou ^ ( ~ -(1 - 3 sin^^i,) - 2 sin^^as - 1 , (26) 

\msj 2 

to second order in the ratio e/m<,. The correlation (see Eqs. (24) and (26)) between the 
solar and atmospheric mixing parameters (sin^6'i2 and sin^6'23) and the active-sterile mixing 
parameter (sin^^i4) is shown in Fig. 5, and is the same for both mass orderings. 

Although this model appears to predict ^13 = 0, this result depends on the triplet VEV 
alignments in the scalar sector. In the general case [40, 41], these alignments could be 
modified to (ip) = {v , v ef" , v ef") and {ip') = {v',v'{l + ei),f'(l + €2)), where the deviation 
parameters come from higher dimensional operators or renormalization group effects. Mixing 
angles will then receive corrections of the same order, proportional to ef^2 ^1,2-, so that one 
can accommodate the latest T2K result as well as the allowed range from Ref. [34]. These 
perturbations would also affect the fourth column of Eq. (17), thus altering the results in 
Eqs. (24) and (26). 

We note here that the small mixing with sterile neutrinos will in general modify mixing 
scenarios. The mixing angles of active and sterile neutrinos are of order e/m^, where e is 
any of the entries (M^^^)/^ with / = e,/i, r. Deviations from initial mixing angles 6'i2, 13,23 
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0.04 



0.03 



0.02- 



0.01 




0.53 



FIG. 5: sin^0i4 against sin^0i2 and sin^023) for both the normal and inverted ordering. The dashed 
(black) lines corresponds to the TBM values of sin^0i2 and sin^023) the solid (red) lines indicate 
the 2(T ranges of the parameters and the (red) square is the best-fit point [7, 33]. 



are then of the same order. The example we have discussed here has the particular feature 
that Ue4 = = Ur4 and that the first and third row of U are identical to TBM. This will 
be different in other cases. 

In a more general flavor symmetry context, it is instuctive to note the origin of the 
eigenvector proportional to (0,-1,1,0)^ in Eq. (22). This arises because of a generalized 
fi-T symmetry. Consider an arbitrary 4x4 Majorana mass matrix M^^^. The defining 
matrix P^r for the Z2 corresponding to fi-r symmetry fulfills P^^ = 1 and the invariance 



condition P^.rM^^'^P^ 



M, 



4x4 



SO that and the resulting mass matrix are given by 



P, 



flT 



/l o\ 
10 
10 
\0 1 / 



and M, 



4x4 



/abb d \ 
■ e f g 
■■eg 
\ ■ ■ ■ m J 



(27) 



The eigenvalue e — / of this mass matrix corresponds to an eigenvector proportional to 
(0,-1,1,0)"^. This Z2 invariance, usually present in A4 models [26], arises by spontaneous 
A4 breaking. Regarding the nearly TBM mixing, there is a second Z2 under which M^^^ is 
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invariant. It is defined by the generator 



sol 



/-I 2 2 \ 
2-120 
2 2-10 
\ 3/ 



(28) 



Additional invariance under this Z2 (note that and Psoi commute) requires f = d — e + b 
and d = g. As a result, the eigenvalue 2e — b — d has an eigenvector proportional to 
(—2, 1, 1,0)"^. Furthermore, it holds that f/gs = U^s = Urs- The main features of Eq. (17) 
are therefore explained by the two Z2 symmetries. In addition, the simplicity of the 
model leads to e = —26, but does not modify the mixing properties that arise from the two 
Z2 symmetries. 

Let us denote with S the upper left 3x3 part of P^oi- With a suitably chosen diagonal 
phase matrix T = diag(l, e~*^'^/^, e*^'^/^), under which the charged lepton mass matrix M^Mj 
is invariant, it holds that = (ST)^ = 1. Hence, S and T generate A4, and the charged 
lepton Z3 defined by T arises again by spontaneous A4 breaking. These well known features 
[26] of models are not altered by our modification. 



B. Two sterile neutrinos 



In order to have two sterile neutrinos in the A4 model, one can simply add a second 
sterile singlet z/^j- before, this sterile neutrino is a singlet 1 under A4 and 1 under Z3, 
and carries the [/(1)fn charge F^,^ = 6. With the additional assumption that the sterile 
sector of the mass matrix is diagonal^, the symmetric 5x5 mass matrix is 



2d 



M, 



5x5 



v 



_ _d 
3 3 
2d 
3 



2d 
3 



e 


f \ 


e 


f 


e 


f 












(29) 



Similar statements about the Z2 invariance of properly extended P^^ and Psoi symmetries 
can be made here (see the discussion at the end of the last subsection). Since / and 
arise in analogy to e and in the one sterile neutrino case, we expect that e and /, as 
well as ■ms2 and m^j, are each of similar magnitude, respectively. In analogy to the case 



^ This can be achieved, for example, with an additional discrete symmetry such as Z2 operating only in the 
sterile sector. 
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discussed in the previous subsection, the mass matrix is approximately diagonahzed by 



U 



+ 



/ 


J_ ^ 

V V o 


o\ 




/o 





e 

rris^ 


f 




1 1 1 













e 

rris^ 




_ 1 1 1 





+ 
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n V oe 
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n n n 
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J 


f ^ 




/■^ \ 
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2 Ui^ + 


"lij / 




















2 >vm|^ ^ 






























3e2 




3e/ 












2 











- 


3e/ 




3f 


) 










2mi^ 





(30) 



assuming that a < rris^ ^ that the ratios e/rus^ and f /nis^ are small. The mass eigen- 
values are 



mi = a + , m2 



3e2 



3/2 



m 



-a + c? , m4 = TUs^ + 



S2 



3e2 



m5 



m 



S2 



(31) 

to second order in the ratios e/m^j and f /rris^. 

Once again, the reactor mixing angle retains its TBM value, and the predictions for 
sin^6'i2 and sin^6'23 are 



sin^^ 



12 



sin^^ 




(32) 



in analogy to the 4x4 case. Using the explicit parameterization (5) of the 5x5 mixing 
matrix, the six active-sterile mixing angles can be approximated by 



sin^ft 



m 



sin 9i 



i5 



f 



m 



S2 



1,2,3) 



(33) 



with the additional assumption that U is real. 

Note that without an additional discrete symmetry the sterile sector of Eq. (29) would 
be a 2 X 2 democratic matrix with entries of order m^-^ ^ — 1 eV. In this case the 
active-sterile mixing would be modified by the presence of an additional large 4-5 rotation 
in the overall 5x5 mixing matrix. 
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IV. REALIZATION OF LIGHT STERILE NEUTRINOS IN SEESAW MODELS 



We now discuss how one can accommodate eV-scale sterile neutrinos in seesaw frame- 
works. We start from the simplest type I seesaw scenario and then present an interesting 
extension to this. 



A. Sterile neutrinos from type I seesaw 

In the canonical type I seesaw model, three heavy right-handed neutrinos 
= i^Ri, ^R2, i^Rs) are introduced and the neutrino mass Lagrangian reads 

-Cm = WMdUr + ]fRMRyR + h.c, (34) 

where the Dirac mass matrix Mo is an arbitrary matrix, while Mr is symmetric according 
to the Majorana nature of right-handed neutrinos. In the basis (z/^,, z/^), the neutrino mass 
matrix is a 6 x 6 matrix: 



and if the entries of Mo are all much smaller than the eigenvalues of Mr, light neutrinos 
acquire masses after the heavy right-handed neutrinos are integrated out, viz. 

~ -MoM^^Ml . (36) 

The mixing between heavy and light neutrinos comes from the full diagonalization of M^, 
and is approximately given by 

R - MoM^^ , (37) 

in the basis where Mr is diagonal. The non-unitary mixing matrix relating the light neutrino 
flavor eigenstates and their mass eigenstates is modified to 

V c^(l-]^RB))U , (38) 

where t/ is a unitary matrix satisfying U'^myU* = diag(mi, m2, ms). Thus, in principle the 
type I seesaw mechanism predicts sterile neutrinos mixed with active neutrinos. However, 
with Mo naturally located at the electroweak scale of 10^ GeV and from the requirement of 
sub-eV light neutrino masses, one typically has that Mr ~ 10^^ GeV. The mass and mixing 
parameters are therefore not suited to explain eV-scale light neutrinos. Possible ways out 
of this dilemma are (i) to lower the overall seesaw scale down to Mr ~ eV [42]; or (ii) bring 
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one of the heavy neutrinos and the Dirac mass matrix entries associated with it down to 
the eV scale. This is in the spirit of the z/MSM [43], where however the adjective "hght" 
denotes keV sterile neutrinos, which can act as warm dark matter particles. 

In the first case, one may use for instance Mr ~ 1.5 eV and Md ^ 0.3 eV, for which 
the light neutrino mass scale lies around 0.06 eV and R ~ 0.2 according to Eqs. (36) and 
(37). This is roughly in agreement with the recently reported range of sterile neutrinos 
from the reactor anti-neutrino anomaly. One would expect three sterile neutrinos in this 
case. In this scenario, there is no neutrino-less double beta decay, since both the active and 
sterile neutrinos contribute to the effective mass (mee), and their contributions are exactly 
cancelled [42]. This is a consequence of the zero in the upper left entry of Eq. (35) and 
the fact that all six neutrino masses are below 100 MeV, which is the typical momentum 
exchange in Oz^/9/3. In addition, neither the baryon asymmetry problem nor the dark matter 
problem is solved in this framework, since the sterile neutrinos are too light. It is also 
questionable why M/j should be so tiny compared to the weak scale, since both arise from 
Yukawa couplings to the SM Higgs doublet. 

Consider now the second possibility, namely bringing one right-handed neutrino down to 
the eV scale. In order to simultaneously explain the baryon asymmetry of the Universe and 
the reactor anomaly, a minimal extension of the SM requires at least one light right-handed 
(sterile) neutrino Vg together with two heavy right-handed neutrinos {vri and i'r2)- In this 
case, the Lagrangian for neutrino masses is 

-Cm = WMdi^r + WMsUs + ^I^MrUr + hj^figJ^s + h.c, (39) 

where Md, Mr and Ms are 3 x 2, 2 x 2 and 3x1 mass matrices, respectively. The full 
neutrino mass matrix in the basis {ul, u^, z/^) reads 
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/ 



(40) 



where we work in the basis where Mr is diagonal (this is always possible by performing a 
proper basis transformation). 

We assume that two of the right-handed neutrinos are heavy (e.g. above the TeV scale or 
higher) and /is is much smaller than Mr, i.e. around the eV scale. Then, at energy scales 
much lower than Mr, vri should be decoupled from the full theory, and one obtains a 4 x 4 
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matrix for active and sterile neutrinos, viz. 



,4X4 _ / -MdM],^MI M; 



In the case where fis ^ oiie may apply the seesaw formula once again, and arrive at 
the light neutrino mass matrix 

~ -MdM^^MI - MsHs^M^ . (42) 

The mixing between active and sterile neutrinos is characterized hy R = MsHg^. Compari- 
son with Eq. (2) gives 

sin ^14 ^ {Ms)iifis\ (43) 
together with the mass of the sterile neutrino 

rris ^ lis- (44) 



The scenario under discussion permits neutrino-less double beta decay, since the active and 
sterile neutrino contributions do not cancel with each other. Furthermore, in this scenario 
the baryon asymmetry of the Universe can be explained via the leptogenesis mechanism, 
for which at least two heavy right-handed neutrinos are required. However, it should also 
be noticed that none of the singlet neutrinos have masses in the allowed range of sterile 
neutrino warm dark matter (i.e. around keV), so that one cannot simultaneously solve the 
dark matter puzzle [43]. 

The question arises how to make one of the right-handed neutrinos so light compared 
to the other two. One obvious possibility is to start with a flavor symmetry that predicts 
one of the masses to be zero. Small breaking of the symmetry will then generate a small 
but non-zero mass. This simple idea has been pursued in the context of keV dark matter 
[44, 45], and will work for light sterile neutrinos as well. One possible symmetry that can 
be used is Lg — — [46, 47]. 



B. Minimal extended type I seesaw 

The cases discussed in the last subsection have the feature that the fundamental seesaw 
Lagrangian already contains at least one particle at the desired eV scale. This is somewhat 
contradictory to the seesaw spirit, and one may ask whether it is also possible to generate 
eV-scale sterile neutrinos without the a priori presence of such states. 

Here we consider an interesting extension of the type I seesaw model, in which three 
right-handed neutrinos together with one singlet S are introduced. A similar idea was used 
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in Ref. [48], with a sterile state of mass ~ 10~^ eV introduced in order to explain the solar 
neutrino problem. Here we will show that there is a natural eV-scale sterile neutrino in this 
scenario, without the need of inserting a small mass term for u^. The Lagrangian of this 
scenario is given by 

-Cm = ITZMdVr + WMsUr + ]^MrUr + h.c, (45) 

where Ms is a 1 x 3 matrix. The full 7x7 neutrino mass matrix in the basis (i^l, i^^, S'^) 
reads 



M, 



7x7 



/ Md \ 

Ml Mr MJ 

V Ms / 



(46) 



In the case where Mr ^ Ms > Md, one should first decouple heavy right-handed neutrinos 
using the canonical seesaw formula, and the effective neutrino mass matrix in the basis 
(z/i, S'^) is given by 

Mr^-f..^SL!!°!!«:5i. (47) 



Ms (M^^) Ml MsM^^M', 



Since M5 is larger than M/5 by definition, one can apply the seesaw formula once again and 
obtain 

^ MdM^^MI [MsM^^MI) Ms [M^^Y Ml - MdM^^MI , (48) 

for the active neutrinos, whereas there exists one sterile neutrino with mass 

~ -MsM^^MJ . (49) 

Note that the right- hand- side of Eq. (48) does not vanish since Ms is a vector rather than a 
square matrix; if Ms were a square matrix this would lead to an exact cancellation between 
the two terms of Eq. (48). 

The active-sterile neutrino mixing matrix takes a 4 x 4 form, i.e.. 



where R = MdM]^^M'^ {MsM]^^M'^) \ Essentially, V14 (i.e. sin6'i4) is suppressed by the 
ratio 0{Mn/Ms). 

As a naive numerical example, for M^, ~ 100 GeV, Ms — 500 GeV and 
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Mji ~ 2 X 10^ GeV, one may estimate that rriy ~ 0.05 eV, ~ 1.3 eV together with 
R ^ 0.2. This is in very good agreement with the fitted sterile neutrino parameters from 
Table I. Neutrino-less double beta decay is again allowed because not all of the neutrinos 
are light. 

Note that one of the light neutrinos is massless since the rank of MJ^^^ is three. This 
model is a minimal extension of the type I seesaw in the sense that one needs at least three 
heavy neutrinos to suppress the masses of both active and sterile neutrinos. In other words, 
two heavy right-handed neutrinos give rise to two massive active neutrinos, while the other 
one is responsible for the mass of Vg- 

Unfortunately, if we embed this scenario into a grand unified theory framework it cannot 
be gauge anomaly free since we have only one generation of S. It is also not possible to 
accommodate two eV-scale sterile neutrinos in this scenario, unless the number of heavy 
neutrinos is increased. Apart from these shortcomings, the scenario possesses the following 
features: 

• apart from the electroweak and seesaw scales, one does not artificially insert small 
mass scales for sterile neutrino masses. As in the canonical type I seesaw, one can 
take Ms > Mr, ~ 0(10^ GeV), while Mr can be chosen close to the B — L scale, not 
far from the grand unification scale; 

• it is a minimal extension of the type I seesaw since three heavy right-handed neutrinos 
can lead to at most three massive light neutrinos (the "seesaw-fair-play-rule" [49]), out 
of which two are active and needed to account for the solar and atmospheric neutrino 
mixing. It is more predictive owing to the absence of one active neutrino mass, while 
it still accommodates all the experimental data; 

• there exist heavy right-handed neutrinos that could be responsible for thermal leptoge- 
nesis. Note that, in the setup we considered, right-handed neutrinos would preferably 
decay to the sterile neutrino since their coupling to S is larger than that to active neu- 
trinos. However, this drawback could be circumvented since 5* enters in the one-loop 
self-energy diagram of the decay of right-handed neutrinos, which could compensate 
for this. 

In this section we have exclusively used the type I seesaw and an extension of it. We 
note that type II seesaw does not provide the possibility of light sterile neutrinos, due to 
the absence of fermionic degrees of freedom. Type III seesaw is formally analogous to type 
I, because the neutral components of the fermion triplets play the role of heavy neutrinos. 
However, these components are not gauge singlets. 



24 



V. CONCLUSION 



Motivated by recently reported hints for sterile neutrinos from particle physics, astropar- 
ticle physics and cosmology, we have studied the realization of light eV-scale sterile neutrinos 
in both model-independent and -dependent manners. The phenomenological consequences 
of light sterile neutrinos were discussed and we pointed out that the presence of such light 
sterile neutrinos would significantly change the oscillation probabilities in short-baseline 
experiments as well as the effective mass measured in neutrino-less double beta decay ex- 
periments. It turns out that the sterile contributions to the effective mass could lead to a 
vanishing (mge) even if the three active neutrinos obey an inverted ordering, which is clearly 
different from the predictions of the standard picture of three active neutrinos. Scenarios in 
which the sterile states are heavier than the active ones are safest in what regards cosmo- 
logical mass limits; the other cases turn out to have rather large contributions, and are even 
in danger of being ruled out by current limits on the effective mass. 

We then focused on the possibility of embedding sterile neutrinos into an effective theory 
with the A4 flavor symmetry, and showed the corrections to the TBM pattern that result 
from this. We found that light sterile neutrinos can naturally be accommodated in the 
flavor A4 symmetry, and their mixing with active neutrinos is correlated to the deviation 
from exact TBM. The departure from the initial mixing scheme is a general feature of 
such approaches; we only studied one example in detail. Furthermore, we have shown that 
light sterile neutrinos can be realized in extensions of the seesaw model. Regardless of 
naturalness, the type I seesaw with three right-handed neutrinos could, in principle, provide 
candidates for light sterile neutrinos. Besides the type I seesaw framework, we have presented 
a minimal extended type I seesaw model, in which, without the need of introducing tiny 
Yukawa couplings, the smallness of sterile neutrino masses is ascribed to the existence of 
heavy singlet neutrinos, whereas the mixing between active and sterile neutrinos could still 
be sizable. Note that all of the models we have discussed share the feature that the sterile 
neutrinos are heavier than the active ones. 

As mentioned before, we only concentrate on eV-scale sterile neutrinos in the current 
study. In general, a keV sterile neutrino, which may lead to rich phenomenological conse- 
quences in the early Universe and supernova explosions, could also be accommodated in the 
models discussed in this work. A detailed survey could be useful and will be elaborated on 
in future. 
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Appendix A: Short-baseline neutrino oscillations: 3+2/2+3 vs. 1+3+1 

When neutrinos propagate in vacuum, the probability of finding a neutrino of initial 
flavor a in the flavor state /3 at the baseline L is given by 



* -\EiL 



(A-1) 



5o.p-AY, Re {U:,U^,Up,U;^) sin^ (^^^ + 2 J] Im (f/:,f/.,f/^.f/;,) sin 



where Ajj = — m'j. The current global-fit of neutrino oscillation experiments indicates 
that (at 2a) Ami = ^21 = (7.24 ... 7.99) x 10^^ eV^ and Ami = l^3il = (2.17 ... 2.64) x 
10"'^ eV^ and that the corresponding mixing angles are sin^ 612 = 0.28 . . . 0.35, sin^ 623 = 
0.41... 0.61 and sin^ 6'i3 < 0.031 [33], where we have neglected small differences between 
normal and inverted ordering, which are irrelevant for our discussion. The parameters 
associated with sterile states are given in Table 1. The oscillation probability for anti- 
neutrinos can be simply obtained from Eq. (A-1) by the replacement U ^ U* . 

It should be stressed that for a terrestrial based neutrino oscillation experiment the 
baseline length is usually fixed, which allows one to make reasonable expansions of the 
oscillation probabilities. In a reactor anti-neutrino oscillation experiment, e.g. Double Chooz 
or Daya Bay, the detector is typically located at a baseline of around 1 km, while the average 
energy of neutrinos is ~ 4 MeV. One can thus estimate that L/E ~ 0{W) eV"^ 
indicating that the oscillation term containing A21 can be safely ignored. In addition, the 
oscillation frequency related to the mass-squared difference A^i {A^i) is much smaller than 
the baseline, and thus the (Asj) term generates a fast oscillation. The anti-neutrino 
survival probability then approximates to 

Pv.-.V. ^ l-2|t/e4|'-4|[/e3|%in2(^^ 

= 1-2 sin^ - cos^ ^14 sin^ (2^13) sin^ ( ) , (A-2) 



AE 



for the case of one sterile neutrino (1+3 and 3+1), and 



P^^^^^ ^ I_2|f/,4|2-2|f/e5|'-4|f/e3|%in2(^^j , (A-3) 

for the case of two sterile neutrinos (3+2, 2+3 and 1+3+1). It is clearly seen that the effect of 
the sterile neutrino (s) is merely to reduce the total neutrino flux, which is needed to explain 
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the recently reported reactor anti-neutrino anomaly. Note that if in the Double Chooz 
experiment the neutrino flux is calibrated using the near detector instead of a Monte Carlo 
simulation, the total shift of the neutrino flux is cancelled when comparing the number of 
events in the near and far detectors, thus having no effect on the measurement of 6*13. 

In order to observe the oscillation effects induced by sterile neutrinos at a reactor neutrino 
experiment, a very short baseline, e.g., L ~ (9(10 m), turns out to be very attractive. At 
such a short baseline length, the A31 contributions to the oscillation probability can be 
neglected, and the survival probability reads 

P.,^., ^1-4 \U.,\' sin^ = 1 - sin^ {29,,) sin^ , (A-4) 

for the 3+1/1+3 case, and 

P..^.. ^ l-4|f/,,rsin^ -4|f/.,rsin^ , (A-5) 

for the 2+3, 3+2 or 1+3+1 cases. Note that the oscillation term proportional to sin^ ('^i^) 
is suppressed by a factor \UeAUe^\ , and can hence be neglected. The practical realization of 
a detector so close to a reactor is not straightforward, and will not be discussed here (see 
Ref. [50] for a recent proposal). So far there is no possibility to distinguish 2+3/3+2 from 
1+3+1 scenarios. 

The LSND and MiniBooNE experiments report an excess of electron anti-neutrino events 
in the oscillation channel I?^ — )■ r'e, with a baseline and energy setup L/E ~ C(l) eV-^ In 
analogy to the very short baseline reactor neutrino oscillations, the oscillation term related 
to A31 can be ignored, and the transition probability reads 

- 4|[/.4nt/,4psin^ (%^) , (A-6) 

for the 3+1/1+3 case, and 

P.,.^.. ^ 4|^e4nf/.4psin^(^)+4|t/.5nt/,5rsin2(%^) 

+ me.U,,U,.U,,\ sin sin cos + , (A-7) 

for the five neutrino case with 5 = arg {UlJJ^/JJer,U*^^ ^ ^14 — ^24 — ^is; where we have given 
the phase in the explicit parameterization of Eq. (5). The additional CP violating phase 
in Eq. (A-7) leads to a CP asymmetry in the transition probability between the neutrino 
and anti-neutrino modes, which is of particular interest in view of the null oscillation results 
from the neutrino running of MiniBooNE. In addition, the last term of Eq. (A-7) indicates 
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FIG. 6: The oscillation probabilities with respect to the quantity L/E m. the 3+2/2+3 cases (solid 
lines) and 1+3+1 case (dashed lines). The mixing parameters are the best-fit values from TABLE 
II of Ref. [33] and TABLE II of Ref. [7]. In the upper panel, we set all CP violating phases to be 
zero, whereas in the lower panel 5 = ^ is assumed. 

a distinctive difference between the 3+2/2+3 and 1+3+1 cases due to the presence of A54. 
Explicitly, in the 3+2/2+3 cases one has A54 = IA51I — IA41I whereas in the 1+3+1 case, 
A54 ~ IA51I + IA41I, which changes the transition probability regardless of the choice of CP 
violating phases. 

We illustrate the oscillation probabilities with respect to the ratio L/E in Fig. 6 for the 
3+2/2+3 and 1+3+1 cases. It is clear that the transition probabilities are quite different in 
the hierarchy schemes, whereas the survival probability is not sensitive to the sterile neutrino 
hierarchies since the A54 contributions are suppressed. Furthermore, if non-vanishing CP 
violating phases are included there are visible differences between neutrino and anti-neutrino 
flavor transitions, implying that a detector at very short distance would be an ideal place 
to search for the CP violation related to sterile neutrinos. 
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